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A Nyquist-Like  Test  for  the  Stability  of  Two-Dimensional 
Digital  Filters 

R .DECARLO,  R.  SAEKS,  \\i>  J.  MURRAY 

Abstract- This  paper  constructs  a Nyquist-like  test  for  the  stability  of 
two-dimensional  digital  filters.  The  test  takes  the  form  of  a continuum 
of  classical  one-variable  Nyquist  plots  parameterized  by  the  elements  of 
the  unit  circle  of  the  complex  plane.  Since  the  parameter  space  is  com- 
pact, the  test  can  be  accurately  approximated  by  a finite  number  of 
classical  one-variable  Nyquist  plots  and  is  therefore  readily  imple- 
mented on  a computer. 


Is  I ROl)l  I TIO\ 

A two-dimensional  digital  filter  is  characterized  by  a rational  transfer 
function  in  two  complex  variables 

B{-~)  (1) 
A(zl,z2) 

where  A(zt.  z2)  and  B(zx.  z2)  are  relatively  prime  polynomials  in  zx 
and  z 2.  l or  the  purpose  of  this  paper  we  say  that  the  digital  filter  is 
stable  if  A(zx,  z2)  * 0 for  r 1 1 ^ 1 and  z2 ; < 1 This  structural 
stability  condition  implies  that  the  filter  is  bounded  input  bounded- 
output  stable,  though  as  recently  shown  by  Goodman  fl|  the  condi- 
tion is  actually  slightly  stronger  Huang  showed  that  this  four- 
dimensional  stability  condition  was  actually  equivalent  to  the 
three-dimensional  condition  that  A(zx.  z2)  * 0 for  izji  = 1 and  i z 2 1 < 
1 or  \Z\  1 ^ 1 and  z2  - 0 which  we  use  as  the  basis  of  our  theory. 

The  key  to  the  formulation  of  our  Nyquist-like  theory  is  the  observa- 
tion that  from  an  abstract  analytic  function  point  of  view  the  classical 
one-variable  Nyquist  plot  is  simply  a method  of  determining  whether  or 
not  an  analytic  function  in  one  variablq  has  zeros  in  an  appropriate  re- 
gion by  plotting  the  image  of  the  function  on  the  boundary  of  the  re- 
gion. To  obtain  a Nyquist  theory  in  two  variables  we  therefore  decom- 
pose the  region  of  C2 . in  which  AiZ\ . z2)  is  forbidden  to  have  zeros  by 
Huang’s  theorem  into  the  union  of  a family  of  one-variable  regions  to 
which  the  classical  Nyquist  theorem  applies  More  precisely,  for  real 
a.  0 ^ a < 2n.  we  define  the  disk  f)Q  in  C*2  by 

= z2 ) l*2»<l}  (2) 

and  we  define  the  disk  D by 

n = {ir , . 0)  in  1 * 1 }.  (3) 
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lig.  1 Nyquist  plot  of  .4(z,  . 0)  for  I sample  1 . 


fig.  2.  Nyquist  plots  for  A(eta . z 2)  for  I xample  I 


Now  . Huang's  theorem  muv  be  restated  as  ‘the  digital  filter  1^  stable  il 
and  only  if  A(zx,  z2)  Iijs  no  zeros  in  the  disk'/)  and  l)  , <*  - a,  2n  " 
Observing  that  each  disk  is  fixed  ir.  one  of  its  variables,  t*  p<  !*  nomial 
A(z  1 , z2)  (restricted  to  any  of  the  above  defined  disks)  1 11  mlytii 
function  of  one  variable  and  hei  e the  classical  N ,uist  test  can  be 
used  to  check  for  zeros  within  the  disk  In  particular.  t(z..  . . ) has 
zeros  in  the  disk  DQ  if  and  only  it  the  Nyquist  pl-  t the  one -variable 
function  Aie,Q,  z2)  docs  not  equal  o:  encircle  zero.  Similarly.  l(Zj. 
z2)  has  no  zeros  in  the  disk  /)  it  and  only  if  the  Nyquist  pint  tor  the 
one-variable  function  A(z  \ . Ol  does  not  equal  or  encircle  zero  Combin- 
ing these  observations  we  obtain  the  following  stability  theorem. 

Theorem 

A digital  filter  characterized  by  the  two-variable  transfer  function 
B(z  1 , z2  ) 

~A(zlVz2) 

where  A(zx . z2)  and  B(zx.  z2)  are  relatively  prime  polynomials  in  two 
variables,  is  stable  (structurally  stable)  if  and  only  if  the  Nyquist  plots 
for  the  family  of  one-variable  functions 

A(c,Q.z2)  0 < a < 2n 

and 

4(z,.0) 

do  not  equal  or  encircle  zero 

Although  the  theorem  formally  implies  that  one  check  a continuum 
of  Ny  quist  plots  parameterized  by  the  complex  numbers  ol  magnitude 
one.  in  fact,  since  this  set  of  numbers  is  compact,  one  can  obtain  a test 
with  arbitrarily  good  resolution  using  only  a finite  number  of  plots.  In- 
deed. in  a somewhat  different  context,  the  authors  have  shown  that  a 
similar  continuum  of  Ny  quist  plots  can  actually  be  reduced  to  a single 
plot  without  inducing  any  error  into  the  stability  test  |2|  The  follow- 
ing examples  are  based  on  a finite  approximation  to  the  continuum  of 
plots  required  by  the  theorem. 

Example  l : Let  the  transfer  function  of  a digital  filter  be 


//(z,.r2> 


1 

1 + 0.25z,  + 0.25z2 


Btzx.z2) 

Aizx,z2) 


(4) 


Step  I Draw  the  Nyquist  plot  for  Aizx . 0).  This  curve,  shown  in 
l ig.  1.  does  not  encircle  zero.  So  we  proceed  to  the  next  step  as  out- 
lined in  the  theorem. 

Step  2 Now  consider  the  family  of  Nyquist  plots  for  the  functions 
A(c,Q.  z2);  0 < a < 2n.  This  family  of  curves  does  not  encircle  “0”  as 
indicated  in  l ig.  2.  Thus  the  filter  is  stable. 

Example  2.  Now  consider  the  filter  whose  transfer  function  is 

1 Biz  j , z2 ) 

//(z,.z2)=  * * (5) 

1 + 0.5Z|  +0.5z2  + 1 2z,  z2  A(2i%22) 
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Fig.  4.  Nyquist  plots  for  A(eia,  z2  ) for  Example  2. 


Step  1 Consider  4(2! , 0).  This  Nyquist  plot  is  illustrated  in  Fie  3 

and  docs  not  encircle  zero. 

At  this  point  no  decision  can  be  made  so  proceed  to  Step  2. 

Step  2:  Consider  the  family  of  functions  A(eiQ,  :2);0  < a < 2n 
yquist  plots  for  some  ot  these  functions  are  shown  m f ie.  4 They 
indicate  that  the  filter  is  indeed  unstable. 
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